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ABSTRACT 

Two planets with orbital period ratio approximately 10:1 have been discovered around the star HD 
83443. The inner and more massive planet, HD 834436, has the smallest semi-major axis among all 
currently known exoplanets. Unlike other short period exoplanets, it maintains a substantial orbital 
eccentricity, e\ = 0.079 ± 0.008, in spite of efficient tidal damping. This is a consequence of its secular 
interactions with HD 83443c whose orbital eccentricity e<i = 0.42 ± 0.06. Dissipation, associated with 
tides the star raises in the inner planet, removes energy but not angular momentum from its orbit, while 
secular interactions transfer angular momentum but not energy from the inner to the outer planet's 
orbit. The outward transfer of angular momentum decreases the tidal decay rate of the inner planet's 
orbital eccentricity while increasing that of the outer planet. The alignment of the apsides of the 
planets' orbits is another consequence of tidal and secular interactions. In this state the ratio of their 
orbital eccentricities, ei/e2, depends upon the secular perturbations the planets exert on each other 
and on additional perturbations that enhance the inner planet's precession rate. Tidal and rotational 
distortions of the inner planet along with general relativity provide the most important of these extra 
precessional perturbations, each of which acts to reduce e\je-2- Provided the planets' orbits are coplanar, 
the observed eccentricity ratio uniquely relates sin? and C = (^/^jX-Ri/ Rj) 5 , where the tidal Love 
number, k 2 , and radius, R\, of the inner planet are scaled by their Jovian equivalents. 

Subject headings: celestial mechanics - (stars:) binaries: spectroscopic - (stars:) planetary systems 



1. INTRODUCTION 

Parameters of the H P 83443 planetary system as given 
by Mayor et al. d2000|) are listed in Table [|. The sub- 
stantial eccentricity of the inner planet's orbit, e\ = 
0.079 ± 0.008, is surprising since other exoplanets with 
comparable semi- major axes have circular orbits. In the 
absence of other interactions, tidal dissip ation woul d erod e 
d as dlnet/dt = -1/r with (Goldreich [l963| , Hut |l98l[ ), 
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Adopting values appropriate for Jupiter: radius Ri « 
7 x 10 4 km, m ass m t « 2 x 10 30 g, tidal Love number 
k 2 « 0.5 (Yoder |l995|), an d tidal quality factor, Q ~ 3xl0 5 



(Goldreich & Soter [1966D , yields t - 3 x 10 s yr. This is in 
striking contrast to the old age of the star. Isochrone fit- 
ting suggests an age of about 6.5 Gyr (N. Murray, private 
communication), which is consistent with the absence of 
chromospheric activity and the small projected rotational 
velocity, vsmi < 2 km/ s. 

Secular interactions with HD 83443c are undoubtedly 
responsible for the pers istenc e of HD 834436's orbital ec- 
centricity (Mayor et al. 200C). Results from calcul ations 
by Lasker and Benz displayed in Mayor et al. ( 2000| ) show 
that secular interactions between the inner and outer plan- 
ets induce anti-correlated oscillations in the two planets' 
orbital eccentricities. The anti-correlation results because 
secular interactions transfer angular momentum, J, but 
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not energy, E, and e is related to E and J by 
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Similar phenomena are noted i n the h ierarchical triple star 
systems HP 1 09648 (Jha et al. |2000D , HD 284163 (Griffi n 
& Gunn 1981 ), and Algol (see, e.g., Kiseleva et al. 1998 ). 

We argue that eccentricity oscillations were tidally 
damped a long time ago, and that currently the planets 
maintain apse alignment (§§). In this state the orbital 
eccentricity of the inner planet is proportional to that of 
the outer planet, and both are steady on the secular time- 
scale. On the much longer tidal time-scale, the individual 
eccentricities decay, with the decay rate being faster for 
the inner planet. 

We utilize the properties of this system to constrain pa- 
rameters that are not directly measurable. From the ratio 
of the two eccentricities, we place bounds on the inclina- 
tion angle, i, of the planets' orbit plane and on the radius 
of the inner planet. 

2. CONSEQUENCES OF APSE ALIGNMENT 

2.1. Assumptions 

We study secular interactions in £fD83443's planetary 
system under the assumptions that the planets share 
a common orbital plane and that their apsides remain 
aligned. The former assumption seems reasonable for plan- 
ets formed in a gaseous disk. The latter assumption is 
consistent with the currently observed approximate apse 
alignment (see Table [l]) and is justified by the analysis in 
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Table 1 

Parameters for the two planets of HD 83443. 



Planet 


P(d) 


e 


w(°) 


a (AU) 


m sin i 


HD834436 
HD83443c 


2.9853 ±0.0009 
29.83 ±0.18 


0.079 ±0.008 
0.42 ±0.06 


300 ± 17 
337 ± 10 


0.0376 
0.175 


1.14 
0.53 



a Data taken from Mayor et al. (2000). Planet masses are in unit of Saturn 
mass M Sat « 5.68 x 10 29 g. The host star A/* = O.79M . 



§||. These assumptions allow us to infer parameters of the 
planetary system from the observed eccentricity ratio. 

2.2. Secular Dynamics 

Mutual secular perturbations of the planets' o rbits are 
calculated following a method devised by Gauss (1818) as 
outlined in the appendix. Equation ( |A5| ) describes the re- 
sulting periapse and eccentricity variations. However, as a 
result of inner planet's proximity to the star, it is subject 
to additional precessional perturbations of which the most 
important ari se fro m its related tidal and rotational distor- 
tions (Sterne 1939] ) 3 , and from general relativity (Einstein 
1916). Expressions for these precession rates are, 
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Note that in equation ( |3b| ) the planet's spin rate, Qi, is 
scaled by its mean motion, n±, since we expect that tidal 
dissipation has produced approximate spin-orbit synchro- 
nization. We neglect smaller precessional contributions 
due to the rotational oblateness of the star and to the 
tidal distortion of the star by the planet. 

Numerically, the total rate of extra precessions amounts 
to, 
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of angular momentum between the planets. Therefore, 
their orbital eccentricities are constant on the secular time- 
scale. Moreover, the eccentricity ratio, ei/ea, is a mono- 
tonically decreasing function of A. To match the ob- 
servationally determined value of &ije.i = 0.19 requires 
A w 9.1 x 10~ n sin « s _1 . This value is greater than 
the precession rate due to general relativity and provides 
evidence for the tidal and rotational distortion of the inner 
planet. Because the inner planet is not observed to transit, 
sini < 0.99, which implies C > 0.9. Nearly pole-on orbits, 
sin i < 0.20, are excluded if we accept that Ri < 2Rj, a 
conse rvativ e upper-limit as judged from Fig. 1 of Burrows 
et al. 2000. This implies mi < 5.45M sat . Figure [j] shows 
ei as a function of C for different values of sini with e 2 
fixed at its observed value. 
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where we set Hi — ni and introduce the dimensionless 

number C = (&2/fc 2 j)(i?i/i?j) 5 , w ^ n Rj an d k 2 j as- 
signed values appropriate for Jupiter. Observations de- 
termine fc 2 j = 0.494, close to that of an n = 1 polytrope, 
&2 ~ 0.51, where as Sat urn has a smaller Love number, 
fc 2 = 0.317 (Yoder |1995| ) because it is more centrally con- 
densed as a result of possessing a relatively larger core of 
heavy elements. 

We incorporate the additional precessional t erm s into 
the secular perturbation equations (see eq. [Al]). In 
a state of apse alignment, there is no secular exchange 

3 The symbol A?2 also appears as the apsidal motion constant in stellar literature in which case it is a factor 2 smaller than the Love number 



Fig. 1. — Values of the inner planet eccentricity, e\, as a function 
of C for different values of sini. The ten solid lines are for sini 
varying in equal increments from 0.1 to 1.0. The horizontal band is 
centered on the measured value of e\ and its vertical width spans 
a 1 — a erroij bar in each direction. Relevant parameters are taken 
from Tabic 



2.3. Tidal Dynamics 
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Tidal dissipation in the approximately synchronously 
rotating inner planet decreases its orbital energy without 
significantly changing its orbital angular momentum. It 
follows from equations (0) and (||) that this causes a\ to 
vary as d\nai/dt « -2ef/r. Since secular interactions 
transfer angular momentum but not energy between the 
planets' orbits, ai remains constant. Conservation of total 
angular momentum leads to, 
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where Ji = yl - efmiriiaf. Apse alignment implies 
e i = /( a ii e 2)i from which we deduce that 
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where f Ul 

merically, f ai ~ 8 and f e2 » 1, so currently ei and e2 
are damping on time-scales of 8.6t and 40r, respectively; 
ei damps faster than ei because the inner planet's semi- 
major axis is shrinking. 

3. LINEAR APPROXIMATION 

In §H we assumed that tidal dissipation had delivered the 
planets' orbits into a state of periapse alignment. Here we 
justify this assumption by studying the secular dynamics 
and tidal evolution under the somewhat inaccurate ap- 
proximation of small orbital eccentricities. 

Our starting point is the expansion of the secular dis- 
turbing function to second order in e\ and e 2 . Substi- 
tuting the expanded disturbing function into Lagrange's 
equations for the variations of the e's and the w's, and 
adding a term to represent tidal damping of e\, yields the 
following set of linear equations (MD2000) 

-j^- = iAixh + iAnh - h/r, ^ = iA 2 ih + iAnh, 
dt dt 

. W 

for the complex variable Ii = e.; exp(in7j). The coefficients 
Aij read: 
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with a = 01/02, Bi = b 
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The symbols A and r retain 



usual Laplace coefficient 
their definitions from §0. 

In the absence of tidal damping (r — > 00) , the general so- 
lution of equation (Q) can be expressed as a super-position 
of two eigenvectors with 

J* = [ei* exp(i e 2 * exp(i 2 *)] exp(ig*t), (9) 
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We distinguish the two eigensolutions by * = p and 
* = m according to the plus or minus sign that ap- 
pears before the square root in the expression for g*. 
Since (ei p /e 2p ) (ei m /e 2m ) = A12/A21 « 1.0, one solu- 
tion satisfies ei < ei and the other ei > ei. If we take 
A = 5 x 10 _11 s _1 , the * = m solution reproduces the ob- 
served eccentricity ratio ei m /ei m — 0.19 and apse align- 
ment. 4 The other solution has ei p /e 2p = 5.26 and anti- 
aligned apses. If both eigenvectors have non-zero lengths, 
the orbital eccentricities oscillate out of phase on time- 
scales of tho usands of years (~ I/5*) as depicted in Mayor 
et al. fl200Cp . 



When r is finite, ai and therefore A^ vary with time. 
Since r 3> 1/fl*, the solution for the eigenvectors given 
by equation ( |10| ) remains valid. However, the eigenvector 
components, e«, decay slowly with time as 
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where 7« = —dlne^/dt. In the limit of small eccentric- 
ity, the expressions for the nonlinear dampingrates given 
by equation (^J) reduce to those in equation ( pr|) provided 
f e2 = 1 as is appropriate for the linear problem. Each of 
the four eigenvector components has its individual damp- 
ing rate. Both components of the p eigenvector decay more 
rapidly than those of the m eigenvector. If the current 
state consists mostly the m eigenvector with a small ad- 
mixture of the p eigenvector, both components in the latter 
would decay on time-scale « r, while the dominant com- 
ponents ei m and ei m would decay much more slowly, as 

llm ~ 7t and 72~m ~ 32r - 

It is notable that the asymmetry in damping rates be- 
tween the aligned and the anti-aligned solutions only ap- 
plies if the inner planet has an orbital period < 5 days. For 
longer periods, the additional precession rates (eqs. Q - 
fl3c| ) are unimportant and (ei m /e 2m ) ~ (ei p /e 2p ) ~ 1. 

4. CONCLUSIONS 

Both the survival of HD 834436's orbital eccentricity and 
the locking of its apse to that of HD 83443c are straight- 
forward consequences of secular and tidal dynamics. The 
ratio ei/e 2 provides compelling evidence for the tidal and 
rotational distortion of HD 834436. It constrains unknown 
system parameters by relating C = 
sini. In particular, we find C > 0.9 



{ki/ku){Ri/Rjf to 



This value of A is different from that required in 82.2 because the linear approximation is inaccurate for e2 = 0.42. 
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Apse alignment indicates that a minimum amount of 
tidal dissipation has occurred over the life-time of the sys- 
tem, taken to be N Gyrs. By requiring that this is at least 
3 times longer than r, we deduce that Q < 5 x 10 5 NC sin i. 
It may also be possible to infer a lower limit for Q by in- 
tegrating the orbits backward in time. 

At the present time n\jn 2 = 9.99 ± 0.06, suggestive of 
a 10 : 1 mean motion resonance. This is probably just 



a coincidence. Moreover, no physically important effects 
would be associated with this resonance should it exist. 
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Peale, and Brad Hansen for helpful discussions. P. Goldre- 
ich's research is supported by NSF Grant AST-0098301. 
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APPENDIX 



SECULAR PERTURBATIONS 
Secular changes of the inner planet's periapse angle, w\, and eccentricity, ex, are described by 
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where the subscript 2 denotes effects arising from secular interactions with the outer planet. Similar equations apply to 
the evolution of w 2 and e 2 but the only significant perturbations are those from the inner planet. Expressions for the 
precession rates due to the tidal and rotational distortion of the inner planet and to gene ral relativity are given in the 
main body of the text. We obtain dwx/dt\ 2 by applying a method devised by Gauss (1818) as outlined below. 

Consider a pair of coplanar rings formed by spreading the mass of each planet along its orbit such that the line 
density is inversely proportional to the orbital velocity. We decompose into radial, 1Z, and tangential, T, components 
the gravitational force exerted on a line element of length r\df\ by a line element of length r 2 df 2 , where / denotes true 
anomaly and r distance from the star; 
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with 9 the angle between r*i and r 2 — r% . We remove 9 in favor of </>, the angle between r\ and r 2 (<fi = f 2 + tu 2 — fx — vjx), 
by applying the geometric relations, |ri — r 2 \ sxxx9 = r 2 sxix(j) and |ri — r 2 \ cos 9 — r 2 cos <fi — r\. Then we integrate over 
df 2 to obtain the components of the total force due to planet 2 on the mass of plane t 1 in the line element r\dfi. When 
this result is substituted into Gauss' form of the perturbation equations (Gauss 1818) and integrated over fx, we obtain^ 
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We have verified that when integrating these equations with parameters appropriate to HD 83443's planetary system, we 
obtain res ults th at agree well with those obtained by direct integration of the orbits using the SWIFT package (Levinson 
& Duncan |99l|). 



5 Our derivation is slightly different from that given in MD2000, but the results are identical. 



